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We show that experiments clearly verify the assumptions made by the first-order two-wave cou- 
pling theory for one dimensional lossless unslanted planar volume holographic gratings using the 
beta-value method rather than Kogelnik's K-vector closure method. Apart from the fact that the 
diffraction process is elastic, a much more striking difference between the theories becomes apparent 
particularly in the direction of the diffracted beam in off-Bragg replay. We therefore monitored the 
direction of the diffracted beam as a function of the off-Bragg replay angle in two distinct cases: 
the diffracted beam lies in the plane of incidence and|[b]|the sample surface normal, the grating 



vector and the incoming beam do not form a plane which calls for the vectorial theory and results 
in conical scattering. 
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I. INTRODUCTION 



Theories to analyze diffraction by sinusoidal structures 
making use of coupled waves have been around since the 
thirties for diffraction of light by standing sound waves 
[Tl[2]- In the late sixties experimental holography became 
popular with the invention of lasers leading to a resur- 
gence of coupled wave theories to describe the diffrac- 
tion of light by plane- wave holograms [3]. Interestingly, 
even much earlier an alternative approach using coupled 
modes (dynamical theory of diffraction, DDT) was fa- 
vored for X-rays [1], which was applied later also for 
neutrons [5]. Though Moharam and Gaylord offer an 
exact solution to the diffraction problem - the rigorous 
coupled wave analysis (RCWA) - [B], and demonstrated 
that the DDT is completely equivalent in its rigorous 
variant [7], still approximate theories are mostly used 
for evaluation of experimental data. Among those the 
first-order two-wave coupling theory of Kogelnik turned 
to the favorite one in holography [3]. Here, the second- 
order derivatives of the amplitudes are neglected on the 
basis of the slowly varying envelope approximation. Its 
secret of success seems to be the simplicity of the rel- 
evant equations and the comprehensive treatment of a 
number of useful cases (transmission-reflection gratings, 
slanted gratings, phase- and absorption gratings) for in- 
Bragg and off-Bragg replay. Despite the indisputable 
merits, soon a variant of the theory was published [S] 
that replaces one assumption made by Kogelnik for off- 
Bragg replay - the K-vector closure method (KVCM) - 
by the beta- value method (BVM) . It was clearly pointed 
out by Syms that from a mathematical point of view 
a class of valid first-order theories exists to which both 
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of them belong [S]. The fundamental difference arises 
from neglecting (KVCM) or properly taking into account 
(BVM) the boundary conditions. A comparison of these 
approximate first-order theories with the RCWA led to 
the conclusion 'that the BVM offers definite advantages 
over the KVCM method' [10]. The BVM later was also 
strongly promoted by the books of Solymar |llj and Yeh 
[12^ as well as by setting up a first-order two- wave cou- 
pling theory for anisotropic media |13j . Any of these 
above mentioned publications focuses on the calculation 
of the amplitudes for the diffracted waves. Both meth- 
ods give similar results for the amplitudes, but strongly 
disagree when it comes to the directions of the diffracted 
beam for off-Bragg replay. 

In what follows we will demonstrate that experiments 
focussing on the direction of the first-order diffracted 
beam falsify the KVCM, as presumed previously in the 
case of X-rays P3, and show that the direction of the 
diffracted beam follows the predictions of the BVM. We 
start with a summary of the predictions given by the 
two theories. Next, the results of diffraction experi- 
ments from an unslanted volume transmission grating 
and a comparison of the experiments with the KVCM 
and BVM are performed. We finish with a discussion 
and conclusion. 



II. FIRST-ORDER TWO- WAVE COUPLING 
DIFFRACTION THEORIES: SUMMARY 

The emphasis of this part lies on working out the dif- 
ference between the two theories for the direction of the 
diffracted beam. The issues related to the differences 
on the amplitudes and phases were already addressed in 
detail theoretically and experimentally , see e.g. Refs. 
[TUl [T5] . The conclusion was that they do not differ when 
obeying the Bragg condition, however, the BVM having 
some advantages over the KVCM for a hologram in off- 
Bragg replay. The first-order diffraction efficiency for an 



unslanted phase grating in transmission is 
r]±i = — i/^sinc^yV^ + ^. 



cs 



(1) 



with v = niTrd/Xy^cu/cs, cji = cosO. Here, ni is the 
refractive-index modulation, d the grating thickness, A 
the wavelength in free space, and 9 the angle of incidence 
in the medium. The parameters cs,£, for each theory 
are defined below in the corresponding sections |IIB| and 
|IIC[ respectively. We show wavevector diagrams yielding 
the expected direction of the diffracted beam for the em- 
ployed experimental configurations and a simple analytic 
formula for the (external) excess angle £, i.e., the angle 
exceeding the value of twice the Bragg angle between the 
zero order and first diffraction orders. The notation is as 
follows: wavevectors inside the medium are denoted by 
q^ in free space by fc; angles inside are and 9b] angles 
outside are 8b, 0, A9 and Z, to be defined in the next 
sections. 



A. Experimental geometries 

Experiments were realized for two different geometries 
that are sketched in Fig. [l] 

[a]. Planar geometry: the wavevector k^ of the incident 
beam, grating vector K and sample surface normal 
N are coplanar in the x — z-plane 

[b]. Tilt geometry: fco find K lie in the plane x ~ z 
whereas N is rotated by an angle Z around the 
grating vector. This is equivalent to the experi- 
mental situation depicted in Fig. 3 of Ref. [TB] and 
treated in a rigorous manner in Ref. |17j . 

Note that for an unslanted grating K ■ N = Q applies. 
For the discussion that follows we assume the sample to 
be fixed in the reference frame while the incident beam 
rotates in the x — z plane in contrast to the experimental 
situation shown in Fig. [T] where instead the sample is 
rotated while keeping the incident beam fixed. 



B. KVCM 

Let us define the parameters for Eq. (fTl) as cs = cr 
and £, — {9 — 9B)\K\d/2. For the wavevector qs of the 
diffracted beam the KVCM simply assumes 



qs ^qR± K, 



(2) 



with qji the wavevector of the forward diffracted 
beam. The corresponding wavevector diagrams for the 
diffraction process are shown in Fig. [2] As pointed out 
earlier [15], Eq. ^ comprises a valid choice to obtain a 
1-D scalar wave equation. However, taking Eq. (pi) se- 
riously, the diffraction would be 'inelastic' for off-Bragg 



replay [TUj, i.e., the wavelength of the diffracted beam is 
expected to change upon rotation. For the parameters 
of the experiment described in section [111] the wavelength 
change at ±20° off-Bragg angle would be about ±10%, 
so that the incident green beam color should change to 
either deep blue or yellow, respectively. The expected 
'excess' angle reads: 



e(Ae) 



n(2sineB -sin AG) 



AG 



n^ - sin^ Ae + (2 sin 6^ - sin Ae)2 
29^. (3) 

C. BVM 



1 



are ^ — 
sin 0)2 with 



For the BVM the parameters for Eq. 
{cos9 — cs)(3d/2 and cg = ^1 — (2 sin 9b 
/3 — 2Tm/X. The latter parameter, which is epony- 
mous to the method, is the modulus of the wavevec- 
tors in the medium with average refractive index n and 
sm9B = ±X\K\/{4:TTn). The Floquet-condition in the 
BVM case requires the wavevector of the diffracted beam 
to obey 



qs^qR±K± AqN 



with 



Aq{9) = /3 (c 



v/l-(2si 



sm t/B — sm ( 



(4) 



(5) 



At first let us consider the planar geometry [a] as defined 
in Sec. |II A| By geometric reasoning and applying Snell's 
law at the boundaries we obtain the (external) excess 
angle e as a function of the (external) off-Bragg angle 
Ae as 

e( AG) = arcsin [2 sin Ob - sin AO] + AO - 26^ ■ (6) 

We want to draw the attention to the fact that e is inde- 
pendent of the refractive index of the medium because of 
phase-matching at the boundaries. The wavevector di- 
agrams for Ae = -22,0, +22° are shown in Fig. [S] It 
can be seen that for qs^x > 27r/A the diffracted beam is 
excited but totally reflected in the medium. The corre- 
sponding rotation angle for which total reflection occurs 
is 



Qtb. = arcsin (2 sin 9s + 1) — 9_b • 



(7) 



Now, let us turn to the more general case, i.e., configu- 
ration [[b]] defined in Sec. HA By inspecting Fig. Hit is 
obvious that now the diffracted beam moves out of the 
X — z plane upon off-Bragg replay. The direction of the 
diffracted beam can be analytically expressed as a func- 
tion of the two external angles 9 and Z characterizing 
the incident beam. The three-step-procedure to arrive 
at this is as follows: define the external wavevector fco 



and the plane of incidence in terms of Q and Z, apply 
Snell's law at the entrance boundary and make use of Eq. 
(El) to derive the internal wavevectors Qr^s- Finally we 
employ Snell's law once more for the diffracted beams to 
end up with qR,s- Again - as expected - the direction of 
the diffracted beams does not depend on the refractive 
index of the medium. The result is lengthy and therefore 
we do not give it explicitly here. We also find an angle 
of total reflection Qtr that depends on the tilting angle 
Z or vice versa a cutoff Ztr.{&) 



■JTR 



(9) 



arccos 



2-\/sin0B(sin9B — sin 6) 
cos 9 



(8) 



D. Comparison 

For the planar geometry both theories correctly pre- 
dict that the direction of the diffracted beam lies in the 
plane of incidence, i.e., kR, N, K,ks are coplanar. When 
applying the KVCM we expect that the diffraction angle 
AkR, ks) > 29b for 9 < 9s and Z{kR, ks) < 29b else. 
On the contrary for the BVM there is always an excess 
angle e > for any deviation from the Bragg condition 
as shown in Fig. [5J For the tilt geometry the KVCM 
again predicts that the diffracted beam lies in the x — z- 
plane, as the surface vector normal does not play a role. 
The BVM on the other hand implies that the diffracted 
beam moves out of this plane for the off-Bragg geometry 
as shown in Fig. |4] Thus for the KVCM we arrive at 



a] see Fig. 



the same situation as in planar geometry 
(dashed lines). When employing the BVM instead, how- 
ever, the diffracted beam is out of the x — z-plane upon 
off-Bragg replay, i.e., the spot moves mostly up/down 
on a screen perpendicular to the diffracted beam for in- 
Bragg. The corresponding situation is shown in Fig. [6] 
{Left). The above mentioned differences as predicted by 
the KVCM and BVM are easy to verify in experiments 
and allow - next to the fact that all processes are elastic 
- to discriminate between the methods. 



III. EXPERIMENTAL 

Light optical diffraction experiments were conducted 
on a Si02 nanoparticle-polymer grating ^18^ with a grat- 
ing spacing A = 2Tr/\K\ = 500 nm using an s-polarized 
Ar-ion laser beam at a readout wavelength A = 5f 4.5 nm 
unless noted otherwise. To prepare the beam properly 
(plane wave), we used a beam expansion system followed 
by a diaphragm with a diameter of about 4 mm. The 
effective thickness of the grating at zero tilt angle was 
found to be do = 58.5 ± 0.05/im by fitting Eq. ([!]) to 
the angular dependence of the diffraction efficiency with 
a coupling strength of v{9b) ~ 1.27r/4. A schematic of 
the setup is shown in Fig. [Tl The sample is placed on 
an Eulerian cradle and rotated around a vertical axis for 



an angle 9 with or without a tilt Z. The screen was 
placed at a distance of 265 mm from the sample, so that 
for Bragg incidence the beam hit the screen perpendic- 
ularly. To monitor the position of the beam we took 
photos of the beam spot on the screen using a digital 
camera (Canon EOS Df 0) in auto-exposure mode. Then 
each photo was analyzed as follows: as a starting point 
we found the position of maximum intensity, followed by 
calculating the spot's center of gravity in a region around 
the maximum with about the size of the beam's diame- 
ter (4 mm = 10 pixel). This was defined as the spot's 
position. The size of the open circles in Figs. [6] and [7] 
reflects the spot's size. 



In planar geometry (case [a] 2 = 0) we observed the 



excess angle £(A9) > upon off-Bragg replay as shown 
in Fig. [7J i.e., the spot on the screen moves in the x — z- 
plane and always in the same direction independent of 
the off-Bragg rotation sense. This behavior is in con- 
tradiction to the KVCM. In checking the validity of the 
BVM we employ Eq. ^ and find that the model fits 
the data without any adjustable parameter. To empha- 
size this we performed the experiment also for another 
read-out wavelength A = 476 nm and the model fits as 
well. While we already falsified the KVCM, the latter 
is a another strong indication for the applicability of the 
BVM instead. 

Next we turn to the even more striking case of tilting 
the grating around the grating vector K about an angle 
oi Z Ri 45° . This results in a larger effective thickness d of 
the grating, thus higher angular selectivity and stronger 
effective coupling as well. This feature is utilized as sug- 
gested in Ref. [12] in recent diffraction experiments from 
holographic gratings using cold and very cold neutrons to 
tune the diffraction efficiency for beam splitters or mir- 
rors [2OH22] . First we checked the angle of total reflection. 
The measured value is (17.3±0.05)°, the theoretically ex- 
pected one Qtr{Z — 45) — 9b = 17.0° according to Eq. 
I8J which fits quite well. Next we monitored the position 
of the diffracted beam while rotating by an angle A9 
off-Bragg. Fig. [6] {Right) shows the experimentally ob- 
served position of the diffraction spots on a screen upon 
variation of A9. The graph compiles all photographs 
for a rotational increment of 2° from A9 = —22° ... 8°. 
This graph impressively demonstrates that the direction 
of the diffracted beam is out of the x — z-plane, thus 
being contradictory to the KVCM. To compare the mea- 
sured data with the predictions of the BVM we fitted 
the model to the experimental data by using a single pa- 
rameter - the tilt angle Z - which in principle can be 
also just measured. A4inimizing the least-squares error 
yielded an angle of Z = 41.7°. The fitted positions are 
given by the solid disks, whereas the centers of gravity 
for each diffraction spot are indicated by dashed circles 
together with the photographs of the diffracted spots in 
the background. The agreement between the data and 
the positions predicted by the BVM is excellent. 



IV. DISCUSSION 

In an excellent review Russell [23: compared the 
dynamical diffraction theory,i.e, a two-eigenmodes ap- 
proach, with the two-wave first order coupling theory 
as derived by Kogelnik and came to the conclusion that 
the expressions for the amplitudes are exactly analo- 
gous. Knowing this we gave a correspondence of the 
relevant parameters, i.e., terms for neutrons [21J to fa- 
cilitate the communication between people working in 
light- or neutron optics. However, it was somewhat as- 
tonishing that the results should be exactly the same as 
'In the modal method, the boundary assumes paramount 
importance, acting as a strong discontinuity(. . . ). In 
the coupled-wave picture the boundary has a very sub- 
sidiary role. '123]. Moreover, for the latter a class of dif- 
ferent boundary conditions are permitted ^ which are 
expected to lead to different amplitudes, too. For the 
KVCM the sample surface normal does not appear in 
the derivation at all, whereas for the BVM phase match- 
ing at the boundary is required exactly as for the DDT. 
So why do KVCM and DDT lead to the same expres- 
sions? In short the answer is as follows: in the deriva- 
tion of the dispersion surface (DDT) the quartic equa- 
tion in the magnitude of the permitted wavevectors is 
reduced to a quadratic one. This approximation is valid 
for small deviations from the Bragg condition. The dis- 
persion surface then consists of hyperbolic branches for 
the permitted wavevectors in the grating region. With- 
out this approximation the equation to solve is in fact of 
4th order and such a deviation is expected to occur for far 
off-Bragg replay. For this reason the approximate DDT 
and the KVCM lead to the identical diffraction ampli- 
tude dependencies. However, when considering the full 
dispersion surface we end up with a different expression 
that is very well consistent with the BVM. A detailed 
analysis will be given elsewhere. 

The normalized diffraction efficiency, i.e., rj{6) / 1]{9 b) ^ 
upon off-Bragg replay for the BVM and the KCVM is de- 
picted in Fig. Is] ( Top). The difference between diffraction 
efficiency according to RCWA and the above discussed 
theories, i.e., Koglenik's KVCM, Uchida's BVM and the 
exact DDT are potted in Fig. Is] [Bottom). One can 
clearly notice that the BVM and the exact DDT yield 
very similar results, in particular when considering the 
zeros of the diffraction efficiency, which are determined 
by the boundary conditions. When now considering the 
case of neutrons (x-rays) the corresponding derivations 
given in, e.g., Refs. PS^-^^T always treat only small de- 
viations from the Bragg condition and thus approximate 
the Ewald spheres by planes. In addition, due to the 
fact that the refractive index is very close to unity, i.e. 
1 ± 10~^, the boundary conditions have only negligible 
influence on the dephasing behaviour of the diffraction 
amplitude. In this case any of the theoretical methods 
describes the experimental results equally well. 

We want to conclude the discussion with a few remarks. 
It was noted in literature that the tilt geometry used here 



can be utilized to achieve polarizing effects, e.g., a polar- 
izing beam-splitter |16j as is obvious already by having 
a look on the rigorous theory [T7] . In the course of the 
treatment given above we did not make any analysis of 
the diffracted beam's polarization state while keeping the 
polarization state of the incoming beam perpendicular to 
the X — z-plane. Yet, it is obvious that whenever the di- 
rection of the diffracted beam for this geometry is needed, 
reasoning on the basis of Kogelnik's KVCM fails, see e.g. 
[28]. 

A second remark relates to a cutoff mentioned in Ref. 
|15) . In this context we want to emphasize that the to- 
tal reflection angle discussed above (see Eq. [t]) is not 
this cutoff. The difference is as follows: in our case a 
diffracted beam is excited in the medium but cannot be 
coupled out via a parallel grating slab. In contrast the 
cutoff angle defines the critical value for which a (higher) 
diffraction order cannot be excited any more (evanescent 
wave). 

The results for the direction of the diffracted beam in 
off-Bragg replay might be important for the implementa- 
tion of a novel beam coupling technique [3S] . The latter 
is an improved variant of the grating translation tech- 
nique to determine the amplitudes and possible phases 
of mixed refractive index/absorption gratings [301 - I32] . 

Finally, it is interesting to note that the kinematical 
theory of diffraction - based on the first Born approxima- 
tion and hence neglecting the effect of multiple scatter- 
ing - already yields the correct direction of the diffracted 
beam, see e.g. [26], and in this respect the KVCM was a 
step backwards. We thus recommend using the BVM be- 
cause it is as easy to implement as the KVCM and gives 
overall correct results. 

Let us pose a final question that might have been al- 
ready anticipated by the reader: Why did this - now ob- 
vious - discrepancy between theory and experiment not 
limit the popularity of the corresponding paper? One 
argument could be that people just employ the theory 
to determine the relevant materials parameters. In this 
case typically far off-Bragg replay does not play a role. 
Another could be the fact that Uchida in Ref. [5] fo- 
cusses on attenuation of the modulation along the thick- 
ness rather than on the different choice of wavevector mis- 
match ('no essential difference'). As noted in Ref. 15j, 
the differences between the exact and the approximate 
theory will be recognized only for 'large deviations from 
the Bragg condition, namely, in the sidelobe structure 
for comparatively thin holograms'. As our nanoparticle- 
polymer composites |18j have a large coupling strength 
one can reach considerable diffraction efficiencies already 
at moderate thickness in the range of 10 /im as opposed 
to the extensively investigated photorefractive electroop- 
tic crystals [TTIIT^ . 



V. CONCLUSION 

It is known for decades that the ampHtude of a beam 
diffracted by a volume grating is usually well described 
by the coupled wave theory using the KVCM (Kogelnik). 
The less prevalent BVM (Uchida) was claimed to be su- 
perior for off-Bragg replay but the diffraction efhciency 



does not allow to discriminate easily between the two 
methods [10] . Here, we investigated the direction of the 
diffracted beam which impressively follows the behavior 
derived from the BVM method - even without any free 
fitting parameter - and disagrees with the KVCM. 
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FIG. 1: Schematic of the setup in an off-Bragg replay difTrac- 
tion experiment for planar and tilt geometry. For the notation 
confer to the text. 
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FIG. 2: KVCM: wavevector diagram for three angles of in- 
cidence: Ae = e - Gs = -22,0, +22°. The full (dashed) 
arrows represent the wavevectors inside (outside) the medium. 
The blue shaded sector covers 2Qb, whereas the hatched sec- 
tor gives the 'excess' angle e. Note, that e < for AO < 0. 
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FIG. 3: BVM; wavevector diagram for three angles of inci- 
dence: Ae = e - Bs = -22,0, +22°. The full (dashed) ar- 
rows represent the wavevectors inside (outside) the medium, 
respectively. The blue shaded sector covers 2&b , whereas the 
hatched sector represents the excess angle £ > 0. 




FIG. 4: 
ration 



: Wavevector diagram for the tilt geometric configu- 
[b] for three rotation angles (AO — 0,±17°) where 



(f denote the corresponding wavevectors. The orange curve 
gives the tie-points of the diffracted wavevectors qs for off- 
Bragg replay. Tilt angle Z = 45°. 
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FIG. 5: Excess angle e(Ae) for the BVM (Eq. |6] solid lines) 
and the KVCM (Eq. [s] dashed lines) for two wavelengths: 
351 nm (thin) and 633 nm (thick). 
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FIG. 6: Left: Diffracted beam on a screen for off-Bragg replay 
according to BVM in configuration [[b]] The screen is perpen- 
dicular to the diffracted beam for in-Bragg.7?ig/if: Measured 
position of the diffracted beams on a screen for a sequence of 
rotation angles O — —22 . . . + 8 with an increment of 2° and 
a tilt angle oi Z = 41.7°. Dashed circles are the centers of 
gravity of the spots, the solid disks are the fitted results. 
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FIG. 7: Measured excess angle £(AO) for the experimental 
data at A = 514 nm (solid circle) and A = 476 nm (dashed 
circle). The solid disks are the simulated values according 
to Eq. ||6|. Note that the model is free of any adjustable 
parameter. 
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FIG. 8: Top: Comparison of the normalized first order diffrac- 
tion efficiency upon off-Bragg replay for Uchida's method 
(BVM, red solid) and Kogelnik's method (KVCM, green 
dashed), cf. also Ref. [TD]. Bottom: Difference of the normal- 
ized first order diffraction efficiency to the two-wave rigorous 
coupled wave analysis upon far off-Bragg replay for Uchida's 
BVM (red solid hue), Kogelnik's KVCM/approximate DDT 
(green dashed line), and the exact DDT without approxima- 
tion (blue dotted line). 



